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Corrections to the Thomas-Fermi-type solution of the Gross-Pitaevskii equation are inevitable in 
order to get correctly the frequencies of the low lying modes out of the Bogolyubov equations. These 
corrections are important in the bulk, too, thus the failure of the Thomas-Fermi approximation is 
not confined to the surface. We discuss this effect quantitatively and consider similar phenomena 
of spin fluctuations in Bose-Einstein condensed gases in an optical trap. 



cd 



> 






I 

o 



% 



Several properties of the recently much studied 
trapped Bose-Einstein condensed alkali metal vapors can 
be reasonably well explained in terms of the Gross- 
Pitaevskii equation and of the Hartree-Fock-Bogolyubov- 
Popov equations (for a review see [|l|). These equations 
yield good approximations for the condensate shape and 
size, the chemical potential, the critical temperature and 
for the excitation frequencies. Note that they do not de- 
scribe the damping of the collective excitations, as the 
Bogolyubov approximation neglects terms of the Hamil- 
tonian that are of third and fourth order in the field op- 
erator (after having used the usual Bogolyubov prescrip- 
tion). Although themselves being approximations (even 
within the framework of nonrelativistic many body the- 
ory), they are usually subjects of further approximations 
to get solutions. This is not as if a numerical approach 
were so difficult, but rather to get a deeper understanding 
of the underlying physics and to clarify the significance of 
the involved mechanisms. The quantum hydrodynamical 
description introduced by Stringari |Q| for trapped Bose 
gases has proven to be important in case of large particle 
numbers. 

At large condensates, and for repulsive interaction (i.e., 
for g > 0, as in case of rubidium and sodium vapours) 
one can apply the Thomas- Fermi approximation for the 
solution of the Gross-Pitaevskii equation (P 
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where 



C = -^A + V{r) - /i + 2.g |$o(f)|' + 2g{^^r)^j{r)) (2) 
2m 

In the above equation ^o(f) stands for the 'con- 
densate wavefunction', V(r) denotes the trap poten- 
tial, the average in the last term on then right hand 
side gives the density of the noncondensed atoms and 



/i is the chemical potential which is connected with 
the particle number and the temperature through N — 



Id' 



In the following we shall 



consider such situations only where $o(^) can be taken 
as real. The Thomas-Fermi approximation amounts to 
the omission of the kinetic term --■^A^o{r). Thus one 
gets the equation 
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Certainly, this equation cannot hold when V{f) > /i, as 
the left hand side of Eq.(Bh cannot be negative. Actually 
it is well known that in the regime V{f) '^ fj, which defines 
approximately the surface of the condensate the Thomas- 
Fermi approximation breaks down, which becomes ob- 
vious if one inserts ^o^ir) into the neglected kinetic 
term. For a spherical condensate this surface layer is 

of width (5 « i (^) rxF where ttf = J^ stands 
for the Thomas-Fermi radius of the condensate and ojq 
is the trapping frequency. Note that ttf oc N^'^ and 
6 ex N~^'^^ ex N~'^'^^rTF- Surface corrections to the 
ground state energy have been thoroughly discussed and 
calculated [|], p, Q. From the point of view of the 
chemical potential and the condensate shape corrections 
to the Thomas-Fermi approximation are important only 
near the surface of the condensate. In contrast, from the 
point of view of the excitation frequencies corrections to 
the Thomas-Fermi approximation are inevitable in the 
bulk as well ||§1 , |P . First of all, recall that at low temper- 
atures the lower excitation frequencies the Hartree-Fock 
Bogolyubov-Popov equations 
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are rather well approximated by the solution of the equa- 
tion 
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called hydrodynamic approximation. For cylindrically 
symmetric harmonic potential (which is relevant for the 
actual experimental situation) and by using the Thomas- 
Fermi expression (j^) at zero temperature, Eq.(|5|) has 
been solved analytically 0, (7|> ^'^"^ ^^"-^ resulting uj 
values agree reasonably with low temperature experimen- 
tal data. Note that at zero temperature Eq. (|l|) is uncou- 
pled from Eq.(^. 

When the radial trapping frequency uoj. is smaller than 
the axial one lo^ (disc-shaped condensate), the lowest 
nonzero eigenvalue in oj — uj^. 

Nevertheless, if one inserts the Thomas-Fermi expres- 
sion for <I>o into the Hartree-Fock-Bogolyubov-Popov 
equation (|4|), one obtains completely erroneous results, 
as we now explicitly show. Indeed, Eqs.(0), (0) imply 
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AF,+2.g(<i>^^) F,=E,G, 



(6) 
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where G,- 



Vj and Fj 



At low energies. 



i.e., when huj <^ fi ^ g ($q ^(0)) the first term on the 
left hand side of Eq. (0) is much smaller than the second 
term, thus can be neglected. Therefore, instead of ^ 



2g($D F,^E,G, 



(8) 



may 



be used. Eliminating Gj one gets 

A($;f^)Vj- = E'jFj. Introducing if = ^l'' Fj 



and (D = -^ one finally arrives at 






(9) 



This equation can be solved with the same techniques 
as Eq.(||), i.e., its validity is extended up to the sur- 
face of the condensate (as defined by the vanishing of 
^'^^) and then it is required that the solution does not 
diverge at the surface. Note that this induces the van- 
ishing of the eigenfunction at the border, unlike in case 
of the hydrodynamical equation where the solution re- 
mains there nonzero. The first eigenvalue at zero tem- 
perature is a; = \/u)1 -f 2ijj'^. One may prove the in- 
equality - /d3r(^$^^ A (^$15^^) > -/d3r(^2^($^^)' 
which readily implies that Q'^ > ujI + 2cu'^. Equality holds 
for the eigenfunction ip = -=<I>J^^. Comparing Cj with 
the experimentally confirmed result uj = cOr the obvious 
conclusion is that Eq.(||) is not correct. We emphasize 
that the source of the error is the Thomas-Fermi approx- 
imation. Indeed, in Ref. |g] Eq.(^ has been derived from 
Eqs.(||) by using the full Gross-Pitaevskii equation (|^). 
That derivation has been done for the zero temperature 
case and for low excitation energies. Taking into account 



the next correction to the Thomas-Fermi approximation, 
one can derive the correct equation (m even at finite tem- 
peratures. To this end, let us write $o 



$r 



$ 



(1) 



where $o^ satisfies Eq.(||) and $0 is the correction 
term. We determine the latter perturbatively (assum- 
ing that we are sufRciently far from the surface), i.e., 
we insert the Ansatz above into Eq.(n]) and keep the ki- 
netic term —H /2m,A$Q ^ while neglect —h /2mA<I>Q 

We get $[)^^ = ^ fttF\i ■ Inserting these 

we get instead of Eqs.(||), (0) 



and U'-Q^) 
expresssions into 
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In Eq.(nl|) we apply the same approximation as in case 
of Eq.(W, namely, we keep only the dominant term 
2g {^0^) Fj on the left hand side, thus we arrive at 
Eq.(||) again. Comparing now Eq.([lO|) with Eq.(^) it is 
clearly seen that the additional term on the left hand 
side cannot be neglected. In fact, for low excitation en- 
ergies the wavelength of the excitations is comparable 
to the condensate size, thus the kinetic energy term is 
of the same order of magnitude as the additional term 
on the left hand side of Eq. (O) B , as we show this be- 
low explicitly for T = 0. Let us emphasize that this 
is true even deeply in the bulk, notwithstanding how 
small $Q is compared to ^q^ ■ Note that in case of 
the Gross-Pitaevskii equation the correction term $q 
should be compared to ^q^ , and they become of the 
same order only in the surface region of the conden- 
sate. Unlike this, when adding the Bogolyubov equa- 
tions to get Eq.(nof), the lowest order terms (i.e., terms 
like 2g (^0^) Gj) cancel exactly, and therefore the cor- 



rection term 



A*: 
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netic term i^AG, 



-Gj should be compared to the ki- 

I^AGj. Eqs.(|l^) and (||) lead readily to 
Eq.(||). Thus taking into account the first correction to 
the Thomas-Fcrmi approximation suffices to recover the 
hydrodynamic approximation, while without it one gets 
the wrong result. We want to investigate this question 
quantitatively. Having solved Eq. (bh at T = and using 
Eq.(ph one can calculate the terms on the left hand side 

(disc-shaped 



of Eq.(nO) separately. In case of ujz > uJr 
condensate) we get for the ratio ^AGj 
the expression 
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where e^ = 1 — co'^/uil and ^, ry are 
spheroidal coordinates defined as p 
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the oblate 



^^e^[e + m-ri% z = / 
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t^r\. For the lowest 



excitation energies to is of the same order of magnitude as 
ujr, therefore within the condensate the ratio of the two 
terms on the left hand side of Eq. (|l^) is of order unity. 

As a further demonstration that the correction term 
describes an important bulk effect we exclude the sur- 
face region, hence we consider a situation when Eq.(||) 
holds for ^ < ^1 < ^0 = i/l/e^ — 1, i.e., inside of the 
condensate, while Eq.(||) holds for ^i < ^ < ^o, i-e., near 
the surface. Note that the Thomas-Fermi surface is the 
ellipsoid £, — £,q and ^ = ^i defines the ellipsoidal border- 
ing surface between the regions of validity of Eqs.(||), (||). 
Then the correction to the eigenvalues of (^ comes from 
the bulk by construction. We calculate this correction by 
perturbation theory in first order, i.e., we consider Eq.(|5|) 
as the unperturbed one and the additional term 



i^'^V^'^V 



(13) 



as a perturbation for ^ < ^i . The correction to the eigen- 
values then has the form 



^ I da d^ie+v')^''^','^'^','' (14) 



Ac^^ 



provided that the unperturbed eigenfunction ip is normed 
by 



" di I driie + v')^' = 1 
J-i 



In case of the lowest excitation frequency one gets 



(15) 



Ac^2 
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3eei(e? + l)[(5-4e2) 



-6^(3 - 2e'M + 1)] + 6-^^i^arctanh (^|^ (16) 

Certainly, the validity of this expression is confined 
for small ^i values only (i.e., perturbation deeply in the 
bulk). In this regime we get in case of w^/'^z ^ V^ 



that ^/ujf+Ku? 



6.74wrCi/Co- The linearity in 



^1 means a proportionality to the volume. Note that the 
considered mode is just the Kohn mode whose frequency 
is known to be exactly the trapping frequency ujr. 

The effect due to the correction term is a consequence 
of inhomogenity (i.e., spatial dependence of ^q^) which 
is induced by the presence of the trapping potential. In 
the Thomas-Fermi approximation the effect depends only 
on the external potential and the number of particles. 
Note that in case of a system contained in a box (with 
zero external potential inside) the Thomas-Fermi solu- 
tion is uniform and the correction of the type discussed 
above is zero. In the extreme situation of a single parti- 
cle, however, spatial dependence of the wave function can 



be a consequence of the boundary conditions, as it is in 
case of a particle in a box. Therefore, it is instructive to 
consider such a surface-induced inhomogenity for a large 
number of particles. In case of a spherical box (with 
zero potential inside) we find that the boundary induces 
a correction to the uniform Thomas-Fermi wave function 
which is comparable with it near the surface and decays 
as exp(— fed) (where k = 2^mfi/h) at a distance d from 
the boundary. Thus, in the middle of the condensate the 
correction term is of order exp(— 1/e), which is nonana- 

lytic in the small parameter ^ = -jtji — \ j^^n (^ being 
the radius of the spherical box). It can be shown that one 
gets a similar nonanalytic boundary induced correction 
in case of smooth trapping potentials, too. 

Let us add that most of the above observations also 
apply to the recently studied Bose condensates in opti- 
cal traps where the Zeeman energy is much smaller than 
the interaction energy. Under such circumstances the to- 
tal spin can freely rotate and new phenomena appear S- 
p2| . Among others, for the total spin / = 1 (as in case of 
^Na and ^^Rb) two types of ground states become pos- 
sible. These are called 'polar' and 'ferromagnetic' [gj and 
possess the form $oC where C is a normalized three com- 
ponent spinor. These ground states correspond to the 
spontaneous breaking of C/(l) x S*^ and 5*0(3) symmetry, 
respectively. For simplicity, we consider a zero temper- 
ature situation. The substitute of the Gross-Pitaevskii 
equation is M 



- — l^ + V{r)-^, + b^l{r) 



$o(r) - 



(17) 



where b = cq for the polar ground state and b = cq + C2 
for the ferromagnetic ground state. Note that cq = 
(.go + 252)73 and C2 = (32 - 5o)/3 with 50 = ATrh^aQ/m 
and 52 = 47r?i a2/m, where qq and 02 are the s-wave scat- 
tering lengths in the total spin F = and F = 2 channel, 
respectively. The polar ground state emerges if C2 > 0, 
while for C2 < the ferromagnetic ground state becomes 
possible. Collective modes can be found in both the po- 
lar and the ferromagnetic ground state as solutions of the 
corresponding generalized Bogolyubov equations. They 
are the density modes in both the polar and ferromag- 
netic states and the spin wave mode in the polar state. 
The mathematical structure of the Bogolyubov equations 
is identical with the scalar case discussed above, there- 
fore, all our previous considerations apply. In contrast, 
the equations, which are related to spin and quadrupolar 
spin fluctuations in the ferromagnetic phase have differ- 
ent structures. The modes have been found in the homo- 
geneous system g and are determined in the presence of 
the external potential by the equations 



- — A + Vir)-fi + g2^l{f) 
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and 
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Here SM'i = 2<^Vi 



It is interesting to look at 



for the spin and quadrupolar spin modes, respectively. 

A combination of Eq.(18) with the Gross-Pitaevskii 
equation ( p^ ) shows that the bulk correction to the 
Thomas-Fermi approximation is important in this case, 
too. Taking this correction into account, we get for the 
spin fluctuations (5A'/_ = <&, 



TF„,Q* 
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= LUjSAI- 
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This is similar to (but not identical with) the hydrody- 
namical equation 0. It is interesting to note that $o('^) 
itself is a solution of the eigenvalue equation (fsl) with 
Ej — (cf. Eq.(p7|)) which corresponds to the solu- 
tion of Eq.(pl (5M„ 



i^V) 



with uji 



0. In case 



of quadrupolar spin fluctuations Eq.(h9f) is relevant. In- 
serting Eq.(O) the lowest order term 2|c2| {^o^) does 
not cancel, therefore one would think bulk corrections 
are negligible in this case, along with the kinetic energy 
of the mode in the bulk. For an extremely large con- 
densate this is indeed the case for the lowest excitations. 
Nevertheless, the numerical value of C2 in case of *^Rb 
(which is the only possible candidate for having a fer- 
romagnetic ground state and / = 1 hyperfine spin) is 
— (9±25) X 10^"^ (72 , very small. Note that if C2 is actually 
not negative, no ferromagnetic ground state is possible. 
The higher order terms (the kinetic term, for instance) 

can be estimated by 2mL'^ ~ -tu ' where L is the major 
axis of the ellipsoidal condensate. On the other hand, 
2|c2| ($0^)^ ~ 2|c2|/i/g2. This will be larger than the 
other terms if 



(21) 



Here ujr is the smallest trapping frequency. For ujr — 
20Hz we get htOr/^i = [10* /N)'^/^, thus for TV > 2 x 10^^ 
Eq.(pll) is satisfied. For larger particle numbers bulk cor- 
rections do not play a role for the lowest lying excitations. 
The kinetic term is also negligible, which looks rather 
paradoxical. The solution is that these excitations are 
confined to the vicinity of the Thomas-Fermi surface as 
noticed already by Ho ||]. For smaller particle numbers 
bulk corrections must be taken into account, and when 
the term 2|c2|$o becomes negligible, the spectra of ( [iq ) 
and (|9|) coincide. Similarly to the case of spin fluctua- 
tions (cf.Eq.(p()|)) one can get the equation describing the 
quadrupolar spin waves in the bulk of the condensate as 
follows: 
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the relative magnitude of the terms in the left hand side 
of Eq.([2^). The second term is dominating if Eq.([2l|) 
holds. In any case its value in the lowest energy state 
(i.e. the gap for excitations confined to the bulk) can 
be estimated by using the Ansatz u'^'~^ = -r^^l)^ ■ One 

obtains for the gap -^^ J cPf(^Q^^ . As the derivation 
is independent of the form of the external potential, this 
result also applies in the homogeneous case. Then we get 
back the known result 2 1 C2 1 1 'I'o P for the gap B . Note that 
for a fixed harmonic external potential the gap increases 
with the particle number as TV^'^ (which is proportional 
to the average density and also t o |< i>(^-'^(0) | in this case). 

It is of worth noting that Eq. (|2l] ) is also the condition 
for the applicability of the hydrodynamic approximation 
in case of the spin fluctuations in the polar state. Due 
to the smallness of C2 Eq.(pi|) is a more confining con- 
dition than the similar conditions in the scalar case and 
in the case of density fluctuations in the polar and fer- 
romagnetic states. Namely, the condition for the latter 
ones is similar to (|2|), but C2 should be replaced by cq. 
This latter condition is also relevant to the applicability 
ofEqs. ^) and (H). 
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